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SUMMARY 


This report presents a theory of oblateness perturba- 
tions of the orbits of artificial satellites based on Hansen's 
theory, with modification for adaptation to fast machine 
computation. The theory permits the easy inclusion of any 
gravitational terms and is suitable for the deduction of geo- 
p ysical and geodetic data from orbit observations on arti- 
lcial satellites. The computations can be carried out to 
any desired order compatible with the accuracy of the geo- 
detic parameters. 
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A MODIFIED HANSEN’S THEORY AS APPLIED 
TO THE MOTION OF ARTIFICIAL SATELLITES 


INTRODUCTION 


This work presents a systematic exposition of a theory of artificial sat- 
ellites based on Hansen's treatment, starting from the basic principles of 
the perturbation theory. 


Hansen's ideas influenced the development of the theories of celestial 
bodies for more than a century; but they are, generally speaking, difficult 
to understand - partly because of the unorthodox treatment of the perturba- 
!!® ns , and pa rtl y because of the way they are presented in his original papers. 
The basic idea consists in introducing a fictitious auxiliary satellite describ- 
ing an auxiliary rotating ellipse of constant shape in accordance with Kepler's 
laws. The position of the real satellite is determined by its deviations in 
time and space from the position of this auxiliary satellite. 

The perturbations of the orbit plane, which are small, are separated 
rom the perturbations in the orbit plane, which are rather large The per- 
turbations of the orbit plane are determined by four interdependent param- 
eters, two of which were introduced by Hansen himself. The perturbations 
in the orbit plane are the perturbations in the radius-vector and in the mean 
anomaly; they are determined by means of a single function W, for which a 
ifferential equation of the first order is formed. The special characteris- 
tics of Hansen s method consist in the addition of the angular perturbations 
to the mean anomaly and not to the true longitude, and also in the use of just 
one unction to determine all the perturbations in the orbit plane. 


In this modification of Hansen's theory the development is a numerical 
one designed to take full advantage of the speed of modern computing ma- 
chines. Repetitive operations also are emphasized for the same reason. 

ma y, m several instances the development has been simplified and short- 
ened by the application of vectorial operations in place of the notations of 
spherical astronomy that were popular in Hansen's time. 
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SOME FORMULAS FROM THE VECTORIAL 
THEORY OF PERTURBATIONS 


It can be assumed in the first approximation that a central body con- 
sisting of homogeneous spherical layers is completely isolated and is rot 
i^g a-und an aL fixed in space. The center of the sphere rs chosen as 
the origin of the inertial system of coordinate s y, z. The * and the y 

axes ar 8 e put in the equatorial plane, and the z axis has the positive direc- 
tion of the axis of rotation. Let i, j, k be the system of basic unit vectors 
connected with this coordinate system. If a satellite wit neg lgi e m 
is placed in the gravitational field of this sphere, the differential equation 
of the satellite’s motion has exactly the same form as m the two-body 
problem: 


r 

r = ' : 

(1) 


(position vector 

r = xi + yj + zk 

of satellite); 

^ i + ^ j + 7- k 
dt dt a- 

(velocity); 

d^x . + d^y . + Jii* k 

dt 2 dt 2 dt 2 

(acceleration). 


The mass of the central body, the radius, and the gravitational constant are 
put equal to unity. The orbit of the satellite is issumed to be an ellipse, 
a more general case it is a conic section. 


The following notations will be used: 


Ei 


a 


E 2 



E 


5 " 1 


semimajor axis of ellipse; 

eccentricity; 

argument of perigee; 

longitude of ascending node in equatorial system; 
inclination of orbit plane to equatorial plane; 
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E. = g Q = mean anomaly at epoch; 



mean motion; 


g = g Q + n(t - t Q ) = mean anomaly; 

E = eccentric anomaly; 

f = true anomaly; 

P = Pi+Pi+Pk = unit vector directed from origin to perigee; 

x z 

R=Ri+Rj+Rk = unit vector normal to orbit plane; 

x y z 

Q = R x P = Q x i + Q y j + Q z k ; 

r° = unit vector in direction of r ; 

n° = unit vector normal to r, lying in orbit plane. 

The first six elements E i (1 = 1, Z, 3, 4, 5, 6) represent the constants of 
integration. The complete solution is given by the following system of clas- 
sical equations: 





P a ( cos E - e 

— E- (Qa/r 

r y/T 


) + Q a \J ^ s i n E ; 

- cos E - P a s ^ n E) * 


(4) 

(5) 


"• ,s»«« r,».v:;'.s s“->- 

f = t° F ( - ! 2 3 4 , 5 , 6 ), which then must be determined by means of 

merits E i \ i - i, 

Equations Z through 5. 

In this exposition, sotne classic r esu.« an. pnes^d in venial 

form. Equation 1 admits of two vectortal integ- als. 


r x r 


( 6 ) 


where 


\f a ( 1 - e 2 ) 

is obtained from Equation 1 by cross multiplying by r and integrating 
result. The Laplacian integral 

Rxr + hr°+g= 0, 


the 


(?) 


where 


is obtained from Equation 1 by cross 
result, taking 


g = h eP , 

multiplying by R and integrating the 


( 8 ) 


r _ _L ( T r - r — ^ = h ^ 
r) x 7? " r 2 \, r ° dt ) dt 

into consideration. Forming the cross produ :t of Equation 7 and R, we have 


R x — - = h(r x 

r 6 


R x (R x f) = 


•R x (hi° + g) . 


RR • f - r 


R . R = -R x (hr° + g) • 


or 
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But 


R * r = 0 , 


and 


R * R = 1 . 


Consequently, 


r = R x (h r° + g ) . (9) 

This is the Hamiltonian integral in vectorial form. The vectors c and g, 
taken as constants of integration, are not independent of each other. The 
condition 


C • g = 0 

must be satisfied, and Equations 6 and 7 supply only five independent scalar 
integrals. 

If the form of the central body is not exactly spherical and the distribu- 
tion of densities is less regular, or if some other bodies — such as the sun 
and the moon — influence the motion of the satellite, then such a motion will 
become a disturbed one and the differential Equation 1 must be replaced by 

r = ~ -qr + F . (10) 


The disturbing force F is assumed to possess a force function ft. In other 
words, we assume that F has the form 


F = grad ft 


3ft . 3ft . 3ft , 

i + j + k . 

3 x 3y 3 z 


(ID 


From the values of r and r for any given moment t , we can deduce 
the system of instantaneous elliptic elements (the system of osculating ele- 
ments) and obtain the position of the instantaneous orbit plane (the osculat- 
ing orbit plane) as a plane passing through r and f: 


r - f (t,a,e,aj,i, t' : , g Q ) 

- f ( t ; Ej , E 2 , ' ' 


r = 0( t , a , e , ^ , i , £ , g Q ) 

= E x , E 2 , 

E,) 


( 13 ) 
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and the elements ate functions of time ” y ‘ he OSC " _ 

lating elements are defined, it can be concluded that 

Br . (14) 

Bt " r ’ 


B r 

Bt“ 


(15) 


sr- mj: -ibr of the 

osculating orbit plane coincides with the rad us vector. We 


Br 

Bt 


Br 

B t 


L = 1 


6 

L 

i=l 


Introducing the differential operator 



Br 

dt 

• B Ei 

1 dE i 

Bi* 

, d t 

■ BE* 

dE t 

B 

dt 

K ’ 


(16) 


(i = 1, 2, 3, 


6 .) 


(17] 


(18) 


1=1 


and taking Equations 14 and 15 into account, we obtain 


and 


b r 

d t 


0 


b f 

d t 


F . 


Cons equently , 



br° 
d t 


0 . 


(19) 


(20) 


( 21 ) 
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These two well known equations serve as a basis for the geometrical theory 
of perturbations. For any osculating element E i (i = 1, Z, 3, * • •, 6), 

SE* dE i 

d t d t. 


In particular, 


dc Sc dg Sg 

dt dt * dt dt 


( 22 ) 


From Equation 6, by applying the S/dt operator, it can be deduced that 

R 4 - T + T- R = rxF. (23) 

dt h h v ' 

T aking 

R • R = 1, R * R = 0 , R x R = 0 , R • r = 0 

into account, we deduce from Equation Z3 that 


-1 I 

d t h 


R * r x F ; 


and 


RxR = hRx (rxF) , 

Rx(rxF) = rR*F“ F R • r . 


Finally, 

RxR = h r R • F , 


and, with 


R x (RxR) - - R 


taken into account, 


R = h(r x R) (R • F) . 
This last equation shows that the vector 


(jJ 


hr (R * F) 


(23a) 
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instantaneous radius vector of the satellite. 

We now introduce a system ° f y Aytogl^hat pla^ 

with the osculating orbit P *”'• ”* Th( , intersect ion of the X-axis with the 
and with the Z- axis norm • • t -/he angle between the X- 

celestial sphere is called t e epar ur p • - Eulerian angles 9, a, i 

axis and the line of nodes is designated by a IhreeJ^ % tem . 

determine the position of the XYZ system with respect to the xyz y 

The satellite is moving in the XY plane. 

in order to investigate the motion of the 
XY7 system we must use the formulas giving th-- conn 

”l„t. y r»d the relative velocities and the relative accelerations. If 
written in the form 

r = Xi' + Yj' + Zk' , 

where i'.j'.k' =R are the basic unit vectors in the XYZ system, then 


dr _ dX , + dY , + dZ 

dT " IT dt J dt 


IS 


the relative velocity, and 


d z r 

d^ 


d^X 

dt 2 


d 2 Y 

+ — 5 " J 
dt 2 


d 2 Z 

dt 2 


k' 


is the relative acceleration. 

The basic formulas giving 
relative motions are 


the connection between the absolute and the 


. dr 

r = + U) x r , 

dt 


(24) 


cHr 

dt 2 


+ 2 U) x 


dr 

dT 


+ OJ x (U) x r) + 


<1(J 
d t 


+ F ; (25) 


or 


, taking Equation 23 into account, 


r - 


dr 

"dt ’ 



3G 


f = 0 + ”( r * £) ,R ' F) 


+ R (R • F) 


- — + F 


By writing F in the form 


F = (F) + — k* 
BZ 


where ( F) is the projection of the disturbing force on the XY plane, 


(F) = 1^- i' + — j' 
BX BY J 


we obtain 


dt 2 = r 3 + (F) • («) 

An important conclusion (Reference 1 ) is noted here: The differential equa- 
tion of the motion of the satellite relative to the system rigidly connected 
with the osculating orbit plane has the same form as in the inertial system. 
The scalar equations of the motion are, consequently: 

In the rectangular coordinates, 


— , x _ Bn 

dt 2 + r 3 “ BX - ’ 


d 2 Y + _Y_ 


and in the polar coordinates, 


(F) = r° + 1 |£ R x r ° , 

or r d v 


v = Z ( X, r) 



(28b) 


(28c) 


Thus, Hansen's classic results about the relative motion can be obtained 
in a more direct way without appealing to the equations for the variation of 
constants (Reference 2). 

Because of the form of Equation 28, which is similar to Equation 1, the 
idea of osculation can be extended to the mov ng system of coordinates; and, 
introducing the operator b/dt, we can claim that in the moving system 

^ = 0, U9> 

tU 

i. i£. r. (F) . (3°) 

d t d t 

The quasi-integrals, Equations 6, 7, and 9 retain their form. In particular 
the area "integral" may be written in the form 


d t h 

The vector H is a constant vector in the XY2 system. By substituting 


R • F 


into Equation 28a we have 


W - hr 


and by substituting Equation 31 into 28c, ancther classic equation is easily 
deduced: 


_d i = m 

d t h ov 


The equation for <lg dt is deduced by apply: ng the operator : .It to the 

Eaplacian "integral": 
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R x — — + hr° + g = 0 . 

at 


( 34 ) 


Taking Equations 34, 33, and 28a into consideration, we obtain 

^7 = 17 r ° x R + (7 + h2 ) IT - r °' (35) 

In the inertial system of coordinates, 

* = fr r ° x R + (7 + h2 ) It r ° + w x K ■ ( 35a ) 

This equation can be used for the computation of the special perturbations 
in the elements. If y is the angular distance of the osculating perigee from 
the departure point, 


g - heP - he ( i / cos y + j ' s in y ) . 

By substituting Equation 36 and 


(36) 


r 0 - i ' cos v + j ' sin v 
r°xR = i 1 sin v - 


J COS V 


into Equation 35, two classic formulas 


d ( he 

c os 
dt 

y) 

_ ^ Bn . 

“ + -=r — sm V + 1 

0 r 

(7 


1 bh 

d ( he 

sin 
d t 

Y) 

BQ 

- “ “7 — cos v + 

0 r 

(7 


an 

B v 


are obtained. These, however, will not be used in this exposition. 


Introduction of the Auxiliary Satellite 


In any satellite theory some first approximation to the real orbit and 
the real motion is used as a starting point. It is customary to call such a 
first approximation an "intermediary orbit." 

The choice of an "intermediary" is not unique. In the theory of Laplace 
as well as in the theory of Hansen, it is a rotating ellipse. In the theory of 
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Hill it is a so-called ’’variational curve,” which in the lunar theory is ob- 
tained by neglecting the eccentricities of the satellite and of the sun and by 
neglecting the solar parallax. In the case of the artificial satellite the "var- 
iational inequalities” depend only on the mean argument of the latitude (Ref- 
erence 3). We imagine a fictitious auxiliary satellite describing the inter- 
mediary orbit in accordance with the prescribed law. The choice of the 
intermediary and of the law of motion on it must be made in such a way that 
the difference between the positions of the real and auxiliary satellites is 
small. 

In the Hansen-type theory which is presented here, the intermediary 
orbit is an ellipse of constant shape lying in the osculating orbit plane, with 
a , e 0 , and n - a “ 3 2 fixed. The auxiliary satellite is describing this 

ellipse in accordance with Kepler's law. The ellipse is rotating uniformly 
with respect to the eccentric anomaly of the auxiliary satellite around the 
axis normal to the osculating plane. The directions and absolute values of 
the radius vectors of the real and auxiliary satellites do not coincide; but the 
difference is small, being of the order of perturbations. Let r° be the unit 
vector along the radius vector r of the real satellite at the time t . The 
radius vector r of the auxiliary satellite will hive the same direction at 
some other moment say ^ . Then, 

r°(t) - r°(j). (37) 

The time t and the "pseudo time” 3 differ from each by the order of pertur- 
bations, and their difference is small. The absolute values of the radius 
vectors r and r also differ from each other by the order of perturbations, 
and we can put 


r ( t. ) = ( 1 + v ) r ( 3 , > , (38) 

where v is small and the factor 1 + v defines how the vector f must be 
"stretched out” or "contracted" in order to become equal with r. Using 
the accepted terminology, call ^ the "disturbeu time," and the difference 
£ - t - £ can be understood as the "perturbation of time.” 

These two simple relations, Equations 37 and 38, serve as a basis for 
Hansen's development. The motion of the auxiliary satellite in the basic 
ellipse is governed by the usual equations famidar from the two-body prob- 
lem. If the true anomaly of the auxiliary satellite is f and the eccentric 
anomaly is E , then 


r c os 


a ( c o s E - t- 0 ) , 



(39) 


r s 1 11 


0 


SIN K , 


(39a) 
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r “ a o (1 “ e o cos E) • (39b) 

E - e Q sin E = g 0 + n 0 * = g 0 + n Q ( t - t 0 ) + <ri Q Sy . (40) 

The area "integral" for the auxiliary satellite takes the form 

o df i : 


d ? 


h > hr 


\/ a n ( 1 - e 2 ) 


It is assumed that the auxiliary ellipse is rotating uniformly with respect to 
E and that the position of its perigee is determined by the equation 

7T = 77 fl + y AE , 


where 


AE = E - E 0 . 

The "angular speed" of rotation y cannot be taken arbitrarily, but must be 
determined in such a way that the development of the perturbations in the 
mean anomaly n 0 8 y does not contain any secular term. The polar angle of 
the auxiliary satellite at the moment 3. with respect to the X-axis is 
f + 7° + ^ AE ' Let the P olar angle of the real satellite at the moment t be 
designated by v. The condition (Equation 37) of equality of the unit vectors 
evidently can be replaced by the condition of the equality of these two polar 
angles: 


V = f + rr 0 + q AE . ( 41 ) 

The node 9 and the argument of the departure point contain two types 
of terms: the secular and the periodic terms. Consequently, the constants 
0 ’ a o * a > and V ma y be determined in such a way that the expressions 

2N = a o + (J o ~ a ~ d ~ 2a AE , (42) 

2 K = cr fl - 0 O - CT- + 0 + 2 r, AE ( 43 ) 

contain neither constant nor secular parts. Such a determination, together 
with the proper determination of y, which is discussed later, leads to the 
development of x, y, z containing periodic terms only. Now 


7 o “ ( a - 77 ) AE - N - K , 


(44) 


/■ i 

0 


- (a + 7;) AE - IS + K . 


(45) 
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For the argument of latitude of the real satellite, 

v - a - f + (tt q - a 0 ) + Uj + a - 7]) AE + N + K . 

The expressions 


(a) 

= a 0 - (a - rj) AE , 

(47) 

(0) 

= 8 q - (a + 77 ) AE , 

(48) 

(«) 

= K - a o ) + + a - V) AE - 

(49) 


containing only the constant and secular parts, are called the mean values 
of the corresponding elements. 


Disturbing Function and its Derivatives 


The disturbing function is defined as the negative of the difference be- 
tween the gravitational potential and the potential of a spherical earth of the 
same mass. In the present theory, ft is taken in the form of an expansion 
in zonal harmonics: 


a 


where 


k k 

— (I - 34> 2 ) + (3>A - 5-^ 3 ) + -J (3 - 30^ z + 35>A 4 ) + ' ' * 

r 3 r 4 r 

n 2 + n 3 + a 4 + - • , 

ip = sin i sin (v — cr) 


(50) 


(51) 


is the sine of the latitude. Taking 

v = f + + V 


(52) 


and 


cr = cr — (a — T)) AE - N — K , 


(53) 


(w) = (77 0 - cr 0 ) + (-y + a - 17 ) AE 


(54) 


into account, we obtain 
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Z - s i n l sin 

[f + ( - ) 

f N + K] . 

( 55 ) 

The expressions for sm N, and K have the form of trigonometric series 

E and (w). It is, however, inconvenient to keep these series in the argu- 
ment in i//, B 

Introducing the parameters 




\ . i 

^ i - sin — cos N , 

A 3 

1 

f-ns j sm K , 


^2 ~ sin sin N , 


<os ~ cos K , 

(56) 


2.2 
- 0 + K 

3 4 

1 , 

(57) 

yields a different form for 0; 




'P = 2(\j\ 4 - \ 2 \ 3 ) sin [f + («)] 

+ 2 (k 2 

4 + ^ i ^ 3 ) ois [f i ( 

")] , (58) 

which leaves only f + («) in the arguments. The parameters X 
were introduced by Hansen. The idea of introducing also and' 

natural one. In the first approximation, ’ 

and a. g 
v, is a 

\ - ■ 1 ° 

- sin 

1 2 ’ 

; \3 ~ 

o , 

(59) 

= ° , 

^4 = « 

1 0 

o s 

2 

(60) 


The next approximations add only small periodic oscillations about the 
60 e only a sHg S htly. ^ V&1UeS ^ fr0m Ec * uations *9 and 

Thus, the introduction of A. parameters permits an easy separation of 
secular motion of the orbit plane from the periodic oscillation of this 

stool To b th S POSi ~° n - mean P ° Siti0n ° f thG P lane is under- 

stood to be the position affected only by the secular motion of the node ) 

about C Te P mea ntS ^ ^ ™ atrix presenting these small oscillations 

. , f ^ n P ° S1 lon of the °rbit plane are polynomials in . The price 

in E q uat ion” 5 1 station of the rotation matrix is ^he condition 

quation 57, that only three parameters, in fact, are independent The 

re ation (Equation 57) is, however, not a great inconvenience and it can be 

ZL three d d^ aCCUraCy ° f the development. If an introduction of 

tory choice woTd IhT 13 nec essary, then probably a satisfac- 
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= 

cos 

N 

i 

tan - , 

(61) 


cos 

K 

2 


It 

1 

1 ^ 

[to 

sin 

N 

i 

t an t , 

(62) 

/V 4 

cos 

K 

2 ’ 


+ ;V 3 

= + tan 

K . 


(63) 


The elements of the matrix of rotation ol the orbit plane about its mean 
position are the rational fractions in p, q, s. The \ parameters are anal- 
ogous to the parameters of Euler in the theory of the rotation of a 
rigid body (Reference 4), and the identity of p, q, s with the components 
of Gibbs' rotation vector (Reference 6) is also easily recognizable. The 
form of v which is convenient to use in developing the disturbing function is 


where 


0 = 2 [(^4 - t + {k 2 k 4 + 


(64) 


/ 3 


s [f + (w)] = 4 ( 1 + n/T - e 2 ) cos [E + (w)] 


+ y (1 - yj~l - e Q 2 ) cos [E - (^)l 


- e Q cos (a;) ; 


(65) 


•m = — sin [f + (u)/j - 0 

a o 1 


)] = -i- (1 + \[\ - e 2 ) sin [E + Ml 


- — ( 1 - J 1 -e 2 ) sin [E “ («) ] 
2 0 


- e Q sin (o>) 


( 66 ) 


For a Q / r 


the classic formula (Reference 6) is: 




(t 


+ 3 cos E - /3 COS 2E + * • 



(G 
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It is necessary to point out that the eccentric anomaly E plays two 
roles. In the motion of the auxiliary satellite in its ellipse, E has the usual 
geometrical meaning; but in the expressions for perturbations it is the in- 
dependent variable replacing the time. These two types of E must be dis- 
tinguished from each other because the partial derivative Bfi/BE, which ap- 
pears in this theory, is taken with respect to the "elliptic" E. The best 
way to keep these two eccentric anomalies separated from each other is to 
use a temporary notation F for the elliptic E in the development of the 
disturbing function. 

The disturbing function can be written in the form 


n ■ = (i + v)'° (i - 3 *a 2 ) 


k a 4 

+ -f (1 + O' 4 <3i/r - 5 1 // 2 ) 
a 4 

n 1 


(1 + O ' 5 (3 - 30 i /> 2 + 350 4 ) — . ( 67 ) 


The elliptic E is replaced by F, and Cl by 
k„ 


Cl' 


(1 + v) 


(1 - 3 i />* 2 ) 


k 3 / a \ 4 

+ — (1 + W' 4 — (3</r* - 5^’ 2 ) 

a n \ P / 


+ — r (1 + O' 5 


(3 - 30'/»* 2 + 3 * 4 ) + 


(68) 


where 


r 


(^1^4 ~ ™* + (^ 2^4 + ^1^3) ^* 


l - — COS [qp + (w)] = -5- (1 + yj 1 _ 


e ^ ) cos [F + (co) ] 


( 69 ) 


o" ( 1 \f 1 - e~^) cos [F “ ( co ) ] - 


e n cos (oj) , 


( 70 ) 
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- sin [<p + (a))] =4(1+ \/l - e^) sm [F + («)] 

a o 

5- ( 1 — J\ - e 2 sin [F ~ (w)] - e sin (w) , (71) 

2 o 

and 

a ° = 2 fi- + 0 cos F + /3 2 cos 2 F + • * * ) . (72) 

* V2 

No replacement of E by F is done in 1 + v, > ^3 ’ ^4 - Conse- 

quently, the introduction of the "temporary" eccentric anomaly F and 0 
instead of 0 permits us to distinguish between two types of E and, if 
necessary, to keep track of changes that occur in fi (or Q*) from one iter- 
ation to another. Each iteration leads to a development of the form 

0* = ECcos (iE + 2j.<*> + *F) + SS sin U E + (2i + D w + * F] 


and 


n = Sc cos (iE + 2}«) + SS sin [-i E + ( 2* + 1) o>] . 

The last equation will not be used. In this exposition there appears the par- 
tial derivative Bfi/BE, which in fact must be understood as 

BO — Bfl (73) 

BE BF 

The "bar" operation means in Hansen's notatioii. the replacement of F by E 
again. Each iteration step leads to an expression of the form 

— = S S sin UE + 2 ->) + SC cos [iE + (2 j + 1) «>] ■ 

BE 


Taking 


3Q _ 3fl 3 r + 3fl 3 f_ 

c F 3 r 3E 3 f 3E 


3 f 


a 0 v 1 e o 



we obtain 


an 

B v 


r an e o sin e 


an 




The development of the derivative 


r (Bfi/Br ) presents no difficult 


le s : 


( 74 ) 


BD _ 3^2 / a n \ 3 

a r ~ (1 + ^)‘ 3 j (i - 30 2 ) 


4k 


f (i + 


*)' 


(30 - 50 3 ) 


5k 


4 c / a 

r (1 + i^)' 5 


*)' 


(3 - 30 4> 2 + 3 50 4 ) + 


= - 3n 2 - 4n 3 - 5 n 4 - • 

= - 30* - 40* - 50* . 


(75) 


Perturbations in the Orbit Plane 

vector ^ 6 Pe + rtUrba * ions in the orbit Plane are the perturbations of the radius 

g ange s variation of elements) are combined into a single angle n 8a 

.hT™ f h f„ e ry PertUrba * i<>nS ° f fhe standpoint^of 
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Differentiating 


give s 


f + r Q + ^ 


dv d£ dj + ii 

dt d? dt dt 


Then, eliminating dv/dt and dt/d, by means of 


d v _ 1 

dt " hr 2 


and 


d£ 

<4 


h o r 


we have 


di = ^ (xY 

dt h \ r j 


-2 


dE 


M h o r dt ’ 


(76) 


or 


, taking into consideration 


r _ 1 

T i + v 


we 


obtain from Equation 76, 


ii = = 

dt h(l + V ) 2 n Q \fl - e 0 2 


r \ 2 dE 
dt 


(77) 


Equation 76 can be written in the form 


d$* : 


dt 


= W + 


> 0 . v' 


h u + v ) 2 „ y i - 


dE 
dt ’ 


(78) 


where 
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W 


1 



1 

1 + v ’ 


(79) 


The first and the third terms in Equation 78 are of the first order, and the 
second term is of the second order with respect to the disturbing forces. 
In order to obtain the main part of & we must find W ; that is the next 
problem. 


The equation of the orbit, 

~ a ( 1 ~ e 2 ) 

1 + e cos f * 


can be written in the form 


1 

h ( 1 + v) 


hr + r * g ; 


and, as a consequence, Equation 79 becomes 


W 


1 - 


2 h Q h r 


2 h o r 


(80) 


By taking 


h 2 
o 


(1 


e o 2 > 


1 


into account, the last equation can be rewritten in the form 


1 - 


2 f 

h o 


(1 - 


% 2) 


+ 2 


h o a o 


(1 - 




(81) 


This is a classic equation written in terms of vectors. In forming the dif- 
ferential equation determining W, it would be preferable to keep the per- 
turbations separated from the elliptic motion of the auxiliary satellite. For 
this reason^it will be more convenient to introduce another function W 
instead of W , by replacing r with p . The vector p is defined as a function 
of the temporary eccentric anomaly F which replaces E and which is con- 
sidered as a temporary constant; or p may also be considered as a function 
of the true anomaly cp , which temporarily replaces f. After the integration 
is performed, F is replaced by E again. Instead of Equations 39 through 
39b, we then have 


P cos cp = a 0 (cos F - e 0 ) , 

p sin Cp - a n V 1 - e 2 
0 0 


(82) 


sin F , 


(82a) 
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P = a 0 (1 " e o cos F) . 


p - i ' p c o s ( v + 77 q + \* - J E ) + J 

In view of the foregoing remarks, W takes -he form 


(82b) 

sin (? + ^ + v-E) . (82c) 


, h o + , h P 

V, = - 1 - IT 2 


2 P ’ g 


(83) 


Substituting Equations 8 
formula: 


W ' 1 h h o a 0 


Eliminating •; in favor of F by means 
obtain another formula: 


h 0 a 0 ( 1 - e o 2 ) h Q a 0 ( 1 - e Q ) 

2c and 36 into Equation 83 results in a classic 

1 + e cos (f + * ’i 11 ^ y ) 


(84) 


1 - e , 


of Equations 82 through 82b, we 


~ + T cos F + ^ iinF, 


(85) 


where 


, _ + o — \ l ~ ee o COi 

- 1 h 2 h fl ! _ e 2 L 


" 0 


^ - x) 


( 86 ) 


+ *.j — E - X') ~ e ( 


1 - e. 


;in (7? 0 + X) 


V = 


(86a) 


(86b) 




It follows from Equations 86 and 86a that 


e o T + = = " 1 h + 2 h • 


(87) 


The formulas 85 and 87 are used in the following exposition in the problem 
a ut e r m i nine the constants of integration and determining the per 


of determining the c 
tions in h/h 0 
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It can easily be seen that W is changing under the influence of the per- 
turbations only. The "elliptic" motion is present in Equation 85 in the form 

of F and is considered as "frozen" for the time being. This leads to the 
conclusion that 


SW dW 

dt dt ‘ (88) 

The perturbations in p are present only through rf + y AE in Equation 82c 
in the argument. The only influence of the perturbative force on p will 
consist of the rotation of p around the Z-axis with the angular velocity 
•y(dE/dt). Consequently, 


s, dE 

dt " (K xp ) q ( 89 ) 

where p is completely uninfluenced by perturbations; then 

a p 

7T “ 0 • (90) 

There is obtained from Equation 83, by applying the operator i/dt and 
taking Equations 88, 89, 90, 35, and 33 into consideration. 


dW 
d t 




o 


P • r° 

a 0 (1 - 


~ P 

e o 2)J 


+ 2h 


o 



r 


3fl 

B r 


+ 


2g • R x p 

h o a o (1 " e 0 2 ) 


M 


dE 

dt • 


Designating by p the unit vector in the direction of p , we have 


(91) 


P = P ■ P° . 

It follows from Equations 82 through 82c that 


±P_ _ ^P_ - d p° d? 

dF " dF H P ~lf IF ' (92) 


dp _ 

W - a o e o sin F 


d'p 

dF 



(92a) 


where 
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Evidently 


f*£l = r x p 

do 


It follows from 92 to 92b that 


dp 

dF 


r = P° a o e o sin F + 


(R x?) ^-x/i - 


Differentiating Equation 83 with respect to F gives 

‘ 


3W 

3F 


_h e o sin F 

h ° 1 - e o 


+ — 


h o a o (1 - e o } 


P° a o e " s 


„ „n F + (R X P ) • -=?- y 1 

0 h 


or 


2g * R x P _ 

p 3 w / 

h 0 \ 

W + — r- + 1 J e sin F 

1 

h o a o (1 - e o 2) 

j; jf \ 

h / 0 

\A - e Q 


Substituting this last result into Equation 91 gives us 


d W_ _ 3H 
dt. o Bv 


P ' r° _ 


1 + 2 


ill P ' r° - p 

h 2 a n ( 1 - e 2 


+ 2h 






5 BW 

o 3F 


-( 


W + — jj- + 1 j e Q sin K 


3r 


dE 
d t 


From Equation 79, 


l o 1 + v 


h 1 - v 

and substituting into Equation 78 gives us 

,2 


(1 + W) , 


dn 0 *> i W + v 

- n A 

0 ’ >2 


i - ^ 


2 \ o 


r 

dt 


(92b) 


• (93) 


dt 


(94) 
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The last equation represents a generalized formula of Hill (Reference 
7) in the theory of planetary perturbations. Differentiating Kepler’s 
equation 


E - e 0 sm E = g 0 + n 0 (t - t 0 ) + n 0 S 9 . 
with respect to E gives 

dt dn 0 S 3- 

a 0 n ° dE + dE 1 (95) 

and, eliminating dt from Equation 94, 

dn o S 9 _ W + v 2 _ /_F_ dn 0 8 l\ 

dE 1 - v 2 \ a o dE ) 

Solving with respect to dn fl S?/dE gives 
dn o S 9 _ w + v 2 ' J _ l - v 2 

dE ' i + w ‘ a 0 ' i + w ’ /fTT 2 ) ' (96) 

0 

Eliminating dn fl S^/dE from Equation 95, 



d t 
*0 dE 


1 - v 2 


’o 1 + W 


— 1 + 


yj 1 - e 2 


(97) 


Taking Equations 74 and 97 into account, we deduce from Equation 93 that 


dW 

dE " h o 


,/TT 


Bfi _ e o sin E 

Si - e 2 3r 


BE 


2 cos (f - cp) - l + 2 

h 2 

0 


P ' COS ( f - Cp) - 1 


1 - e 2 
0 


r 1 - v 2 / 

l o n o 1 + W 1 + 


/TT 


+ 2h 0 — sin ( f - <p) 


p BW 


r ^ 1 - v 2 / 

‘o n o 1 + W ( ^ 


/ 


' r 


1 — e 


Bn 

Br 


A “ e 2 


o L - 


a Q * BF \ W + 1 + h ) e o sin F 


(98) 
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or 


after expressing the results in terms of E and F, 

s-y 


dW = L ^ 3 i_ + Nr -J—l A + 


dE 


BE 


Br 


/TV 


( 99 ) 


where 


M ( 1 - eO 


h 2 
o 


2 + 2e . cos E + 2 cos (F E) 


(F ~ 2E) - e cos F 


" (" 2e 2 

1 + v L 0 


- 2e Q cos 


E + e 2 cos (F + E) + C: - e Q 2 ) cos (F - E) 


- 2 e Q cos F 


+ I - J e 0 2 + 2e 0 cos E 


- 1 — e 2 cos 2E^ , (100) 

2 o / 


,2 r 


NU - = 77 


2e„ sin E - e sin 


F + e fl sin (F - 2E) 


1 + v 


-2e 0 sin E “ (2 - e fl 2 ) sin (F " E) 


>] * (. 


+ e 2 sin (F + E) | + (e„ sin E - y e Q 2 sin 2E^ , (101) 


1 -JL 1 A + — • * — 

1 + M 30 i/l - e n 2 


and 




In the first approximation, there is 
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M ( 1 - e ^ ) 
o 


3 * 1 2 e o 2 


+ 2e 0 cos E “ — e Q 2 cos 2E 


+ e* c 
0 


os (E + F) + ( 4 - e 2 ) cos (E - F) 


-0 cos (2E - F) - 3 e cos F 


N(l-e 2 ) - + e sin E ~-i-e 2 sin2E + 


2 o 


e z sin (E + F) 


+ (2 - e 2 ) sin(E - F) - e sin(2E - F) - e„ sin F 


and 


S - “ 2 sin F. 

Before integration of Equation 99 is started, -y must be determined in 
such a way that no term of the form A sin F is present in this equation. 
Otherwise, the integration will produce a secular term AE sin F in W and 
a term AE sin E of the mixed type in W. As long as the inclination is not 
near 63.4 degrees, only periodic terms are permitted and small divisors 
will not make the process inoperable. 


Perturbations of the Orbit Plane 


It already has been pointed out that the motion of the orbit plane can be 
decomposed into the secular motion of the mean orbit plane and the small 
oscillations about this mean position of the orbit plane. In order to deter- 
mine the secular motion as well as the small oscillations of the orbit plane, 
the differential equations for A parameters must be formed. The classic 
equations for the variation of constants will be used: 


■ de 



dt 

h r sin (v - a) , 

(102) 

dcr 

dO 


d t 

. cos i , 

dt 

(103) 

di 

BO 


d7 = 

h r co s ( v - a) . 

(104) 
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Differentiating Equations 42 and 43 gives us 


dN 

2 dT r 

dcT 

” ^dt 

d O 
dT 

dE 

" 2 ' dT ■ 

(105) 

, i*K 

~ dt 

da 
d t 

d(* 

+ — 
dt 

+ 2,^ . 
1 dt 

(106) 


Substituting Equations 102 and 103 into the twc last equations, we have 


dN 
d t 

dE 
d t 

T- 

3H i 

^ cot Y 

Sill ( v - o ) , 

(107) 

dK 

dt ” + 7j 

dE 

dt + 

7 hr 

Bfi i 

BZ tan T 

sin ( v — a ) . 

(108) 


Differentiating Equations 56 for A parameters and substituting Equations 
104, 107, and 108 yields the following equations: 


d\ 


dt 

dt 

dV 3 

dt 

dt 

But in this case 


+ a 4 ^ + “ a n h ( 1 + v) rq- ( + k A •'£ - A, *0 


dt 2 0 


dE X 

"l dt 2 


3Z 


a 0 


a-\, — + -r a n h(] + v) - \ 4 "i) 


+ V A , + 4" a Q h ( 1 + ^ ( + A. 2 ^ + ^ l 


4 dt 2 


- V 4 ^ + | a 0 h (1 + v) (- q 4 + \ 2 *n) 


"3 dt 2 


(109) 

(109a) 

(109b) 

(109c) 


BO _ B«_ 
3Z dip 


and 


d t 
dE 


A . 


( 110 ) 


( 111 ) 


Taking Equations 110 and 111 into consideration, the following final 
equations can be deduced from Equations 109 through 109c: 
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dE 

d \ 2 

~dE 

dX 3 

~dir 



+ 7 ^4 + i 

~ Vk i + J 


h 

a o 

30 

h o * 

y/l - e 2 

■ w cos 


v 0 


h 

a o 

30 

h o ’ 

\J 1 - e 2 

30 cos 


0 


A 

a o 

30 

L 


COS 

h o 

s/l - e 2 

p 


0 


h 

a o 


h o 

A - e 2 
0 

30 cos 




i (- Ag "C - /V 4 tti) A , 


1 (+ k 2 l + m) A , 


1 ( - A | "t + A 2 171 ) A . 


( 112 ) 


(112a) 


( 112 b) 


(H 2 c) 


The values of a and r) must be determined from the second and the 
fourth equations in such a way that no constant terms are present on the 
right-hand sides. The system (Equations 112 through 112c) is solved by the 
method of iterations starting with 


k 


1 


sin 



2 


= ° , 


Each iteration step leads to an improved value of A. and a, 7 j, and the proc 
ess must be continued until the final values of A. , A. \ A a an( j v 
are reached. 3 ’ 4 ’ ’ 


ON THE CONSTANTS OF INTEGRATION, 
DETERMINATION OF THE BASIC HANSEN FUNCTION, 
AND PERTURBATIONS IN THE MEAN ANOMALY 

The constants of integration are the six given elements a Q , e 0 , g , 6 , 
1 o> anc * = n o ~ a o ‘ moment of time exists for which these elements ° 
are osculating. They must be compatible with the observations, and the 
development of the coordinates obtained on the basis of these elements must 
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contain the periodic terms only. The system oi basic elements is determined 
by repetition of the orbit correction; after each orbit correction the develop- 
ment of perturbations must be done anew in order to obtain a better repre- 
sentation. If the two theories are compared, the systems of elements used 
in different theories are similar but not identical. Such a comparison re- 
quires, generally speaking, great care and ingenuity. But in any given theory 
there is a unique system of elements compatible with the observations. 

In the integration connected with determining the perturbations, some 
new additional constants of integration are introduced. They must be deter- 
mined in such a way that no secular terms appear in the development of the 
coordinates. Implicitly, these additional constants are determined as func- 
tions of the basic elements. The analytical form of the dependence of these 
constants on the elements is not evident in the numerical theory. 

In this development there will be series of the following four forms: 

S C cos (-1 E + 2 j.oo) , 13 > 

SC cos ['i E + (2 i + 1)"] , (113a) 

S S sin (o-E + 2 3 . 0 ;) , (113b) 

L S sin [/IE + (2 j + 1) to] . (113c) 

In order to keep the secular terms absent in the development of the co- 
ordinates, no constant of integration is added if the integrated series has 
the form of Equations 113a, 113b, or 113c; but cne is added if the integrated 
series has the form 113. That is clear because series in Equations 113b 
and 113c, being the sine series, do not contain any constant terms. The 
series 113a cannot contain any constant term because only the combinations 
±co 9 ±3 co, ±5gj, * * * can be present in the argurr ent. 

At each step of the process of iteration, \*e obtain for Equation 99 an 
expression of the form 

= E S sin ( a , E + 2 ^ oj + & F) 


+ 2 C cos UE + ( 2 i + l)co + k F] (114) 


or of the form 
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dW 


“ S S sin (i-E + 2 j.a>) + S C cos [4 E + (2 i + 1 ) co] 


S S sin { -i E + 2 ico) + 2 C cos [-i, E + (2 > + 1) co] 
SC cos ( -i E + 2 yco) + 2 S sin [-t, E + ( 2 -j, + 1 ) o>] 


cos F 


sin F. (11 4a) 


If the series in Equation 114 or 114a are integrated, the additive constant of 
integration must be of the form 


C Q + Cj cos F 

and the series for W has the form 

W = C # + Cj cos F + SC cos('lE + 2 -j, co + 4 F) 

+ S S sin [A, E + ( 2 3 . + 1 ) + &F] . (114b) 

The series for W becomes 


w C 0 + Cj cos E + S C cos (4, E + 2 j,co) + S S sin [A. E + ( 2 3, + 1 ) co] . (114c) 
Equation 96 for n Q 8 ■j, can be put in the form 


dn o S 9 

dE 




1 + W 


r 





1 + W 



(115) 


This form is convenient for the use of the method of iteration because the 
values of v , W 2 and 1 + W can be taken from the previous iteration. 
Hence, Equation 115 takes the form 


dn 0 S 9 

jg - + ^1 cos E) ( 1 - e 0 cos E) 


+ ^ cos + 2 i<y) + S sin [A.E + (2 i + 1) co] . (115a) 

The coefficients C Q and C : must be determined in such a .way that no con- 
stant term and no term of the form K cos E are contained in the last equa- 
tion because such terms, when integrated, will produce terms already con- 
tained in Kepler's equation. Then 
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c 0 - y e o C l + A 00 0 ’ 

- e 0 C 0 + C : + A 10 = 0 , 


or 


2A 00 + e A 10 
2 - e„ 2 


2A 0 o e o + 2A 10 

c - 

and the coefficient A 20 is corrected by an amount 

1 r 

" T 0 1 ‘ 

The final forms of W, W and dn 0 & 7 /dt, after the constants C Q and Cj are 
determined and substituted into Equations 114b, 114c, and 115, remain t e 
same as before: 

w = s c cos UE + 2 ico + h . F) + S S sin Me + (2 i + 1) <*> + &F] . 
w = sc cos(iE + 2 irC*}) + S S sin [iK + (2 i + 1) <A , 


dn ° 8 - = s c cos (iE + 2 j«) + S s sm [A. •: + (2 i + i) ; 
dt 

and at each step of the iteration process we obtain for n fl 8-j a result of the 
form 

n 0 8? = 2 S sin UE + 2 i<w) + S C :os [-CE + (2 i + 1) cu] . 


Determination of the Perturbations in h/h 0 
and in the Radius Vector 

We have obtained for W: 

w = S C cos ( i E + 2 i-u + &F) + 2 S sm [A E + (2 i + 1) w + kf] 
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where 


1 = - 1 , 0 , + 1 . 

This result can be presented in the form 

W = E + T cos F + 'P sin F ; (116) 

where E is obtained as the part in Equation 114b independent of F, and y 
can be obtained by putting F = 0 in the remaining part because 

r = (* - E) f=q . (117) 

The perturbations in h/h 0 can be obtained from 

E+e 0 r = -l--^- + 2h (87) 

h h o 

Putting 


h 

we deduce from the last equation 


+ A , 


A 



T) 



A 


2 


(1 + A) 


(118) 


(119) 


or 


A 




e v ) 
1 ' 



- A" 


•••) • 


( 120 ) 


Also, 




e v ) 
0 A } 


Equation 79 can be put in the form 


( 121 ) 


v - |-(A - W) - ~ v (A + W) . (122) 

This form ^s more convenient for the use of the process of iterations be- 
cause A, W and u in the second term can be taken from the previous 
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iteration. The determination of h/h 0 and 1 + v, consequently, does not require 
any additional integration after the basic function W is determined. 


Determination of the A Parameters 


After the constants a and tj are determined, Equations 112 through 
1 1 Zc give 


sin + 4- (A + B) + 5A.J , 


A , 


5\ 3 , 


(123) 

(124) 

(125) 


The terms 


cos g + 2 ^ B) + • 


s\ 1( sa 2 , 5\ 3 , s;. 4 


( 126 ) 


are obtained by the formal integration of Equations 109 through 109c. They 
do not contain the additive constants of integration, and they have the form 


S\ 1 = SC cos(iE + 2}^) + S S sin [i,E + (2j + l)w] , 
= S S sin ( -i E + 2 ')(<&) + SC cos [a. E + ( 2 q. + 1) , 


+ S C cos [iE t ( 2 -j. + 1 ) w] > 
SC cos ( 't E + 2 y w) + SSsin[-tE+ (2^ + l)c»j] 


S\ 3 - S S sin (-iE + 2 


&k, 


(127) 

(128) 

(129) 

(130) 


Only the series for and A.^ contain the additive constants of integration. 

The determination of these constants can be dene in a more symmetrical 
way if they are written in the form (A + B)/2 and (A - B)/2 . 

Two conditions must be satisfied: (1) The principal term in the latitude 

must have the form 



35 


(Z) In addition, 


Now, Equation 58 is 


sin i Q sin [f + (co) ] ; 


X ’ 2 + k 2 + k 3 + k l 


"1 


(131) 


( 132 ) 


0 - 2(\j\ 4 k 2 k 3 ) sin [f + (<u)] + 2 (k 2 k i + X.^3) cos [f + (*))], ( 58 ) 

and, as a consequence of the first condition, the constant part in 

2(^1 A. 4 - k 2 k 3 ) is equal to sin i Q . Substituting Equations 123 through 126 
into 131 and 13Z, gives us 

* .i. A . (c„. h.. .i„ 

+ const, m 2(Sk } Sk 4 - Sk 2 Sk 3 ) = 0, 

and 

T <a2 * b2 » * (— T + *i» ^)a - (». h.- ,i„ i) B 

2 2 9 n 

+ const, in (SA.J + S\ 2 + SA _3 + 8k ) = 0 , 

or 


where 


A 2 + 2A |cos -j + sin -j-j + (ID = 0 , (133) 

B 2 - 2B ^cos ~ - sin ) + (12) = 0 , (134) 


(11) = cons t , 

in 

(8\j 

+ 8A 4 ) 2 + 

(sa 2 

- SX 3 ) 

( 12 ) = cons t . 

in 

r 

(S\j 

- 8A 4 ) 2 + 

(Sk 2 

+ S V 


Equations^ 133 to 134 are solved by the method of successive approximations 

l / i 0 T 90 . It is preferable to use the p , q , s , system if the satellite is 
a polar one. We simply put 
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x o j. 

p = tan — + dp , 

q = S q , 

s = 8 s; 

and 8p, 8 q, 8 s are deduced by formal integration. 


Decomposition of the Matrix of Rotation 


In his lunar theory, Hansen developed the radius vector, the longitude, 
and the latitude into trigonometric series. In this case, however, the incli- 
nation can be large, and it is preferable to use the matrix of rotation instead 
of the development of the coordinates into senes. Two types of rotation 
matrices appear in the present exposition: 

(1) The matrix representing a rotation about the x-axis for the angle 
a . This matrix has the form 


A l [a] 


"+1 0 0 
o + cos a - sin a 

0 + sin a + cos a_ 


(Z) The matrix representing a rotation about the z-axis: 

+ cos a - si i cl 0 

+ sin a + cos cl 0 

0 0+1 


A 3 [a] 


The coordinates of the real satellite with respect to the system rigidly 
connected to the moving ellipse (the x-axis is directed to the perigee of the 
moving ellipse and the z-axis is normal to tie orbit plane) are: 


(l + v) a 0 ( COS E - e 0 ) , 


( 1 + V) a 0 y 1 - e o 2 


sin E , 


0 . 
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The transformation from the system rigidly connected to the moving 
ellipse to the inertial system requires a triple rotation, and the position 
vector in the inertial system can be presented in the form 


r 


(i + v) r ■ 


a„ (cos E - e Q ) 
a 0 \/l ~ e fl 2 sin E 


r A 3 [9] • Aj[i] ■ A 3 [- 7 t 0 + y AE - cr] 


(135) 


(136) 


We had 


6 = (6) - N + K , 

+ AE “ a = M + N + K , 

and we deduce 

r = A 3 [((9) ] ' A 3 ^ k “ N] • Aj[i] • A 3 [K + N] • A 3 [(o))] . 

Taking 

^1 “ s * n 2" cos ^ f ^3 = cos ~ sin K , 

*2 = sin y sin N , A 4 = COS y cos K 

into account, we have 

r = • [A-jj] • A 3 [ (ccj)] , (137) 

where 


*11 

= * M 2 ■ 

2 

^2 

. 2 , 2 
" *3 + *4 

j 

>* 

to 

II 

+ 2 

(* 3 

*4 - 

*1 * 2 ) , 

*12 

- - 2 (\ 3 


*1 *2^ > 


*22 - 


2 

1 + 

*/- 

, 2 

- * 3 + * 

*1 3 

= + 2 (A.j 

^3 

*2 *4 1 . 


*2 3 = 

- 2 

(*1 

*4 + 

* 2 *3 ) , 



*31 

= + 2 <*3 

*1 + 

A 2 * 4 ^ 

1 






*32 

= + 2(\ 

*1 - 

* 2 *3) 

f 






*33 

= - - 

^ 2 

, 2 

+ * 3 + 

< 

. 
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The numerical computation of the perturbations for a given moment of 
time is done by the method of iterations. Using Kepler’s equation (Equation 
40 ) and the development of the values o ' E and n 0 Sy are obtained by 

the method of iterations, starting with - 0. 


c 


Then the values of ^ij> T and 1 + u arc e 
oordinates x, y, z in the inertial system a -e 


valuated and, finally, the 
evaluated. 


CONCLUDING REMARKS 


The theory described here is a numerical one and permits the full use 
of the large capacity of modern machines. The computation can be carried 
out to any desired order compatible with the accuracy of the basic data. 
Terms are retained or rejected on the basis if their numerical values and 
the decision about the importance of a certain term is made by the machine 
automatically. Hansen’s theory permits the easy inclusion of any number of 
gravitational sources, and the present program can be used without modifi- 
cation for such cases. 

However, the numerical treatment has disadvantages in certain cases. 

If the eccentricity is small, the determination of the motion of the perigee 
is difficult because the eccentricity appears is a divisor. In the case o 
large eccentricity, difficulties arise through the presence of the factor 
1 _ e 2 in the denominator and by the slow convergence of the series or 
a /F ° For these two extreme cases an analytical development would be 
preferable. In this connection we refer to the results obtained by Brouwer 
(Reference 8), Kozai (Reference 9), and Garfinkel (Reference 10). 

In addition, there is the special problem involved in the treatment of 
the critical inclination. An extension of the oresent methods designed for 
the solution of the critical inclination as well as the cases of very small and 
very large eccentricity, is now in progress. 
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